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Optical trapping techniques allow for the formation of bosonic condensates with internal degrees 
of freedom, so-called spinor condensates. Mean-held models of spinor condensates highlight the 
sensitivity of the quantum phases of the system to the relative strength of the two-body interaction 
in the different spin-channels. Such a description captures the regime where these interactions are 
weak. In the opposite and largely unexplored regime of strongly correlated spinor condensates, 
three-body interactions can play an important role through the Ehmov effect, producing possible 
novel phases. Here, we study the three-body spinor problem using the hyperspherical adiabatic 
representation for spin-1, -2 and -3 condensates in the strongly-correlated regime. We characterize 
the Ehmov physics for such systems and the relevant three-body mean-held parameters. We hnd 
that the Ehmov effect can strongly affect the spin dynamics and three-body mean-held contributions 
to the possible quantum phases of the condensate through universal contributions to scattering 
observables. 

PACS numbers: 67.85.Fg, 31.15.xj, 34.50.-s 


I. INTRODUCTION 

Over the last few decades technological advancements 
in trapping have allowed for the creation of a range 
of ultracold dilute gases in the quantum degenerate 
regime. In bosonic systems, quantum statistics lead to 
a global phase coherence and macroscopic occupation of 
the ground state of the system. In spinor condensates 
the atomic spin degrees-of-freedom are energetically ac¬ 
cessible leading to a unique scenario in which both su¬ 
perfluity and magnetism can exist [I|. The first spinor 
condensates were made with ®^Rb in a double magnetic 
trap Q and with ^^Na in an optical trap Q. There has 
been a subsequent explosion of experimental [Ml , and 
theoretical [ihl-l^ studies in the area. Studying the in¬ 
terplay between superfluidity and magnetism has proven 
to be a rich area for probing the many-body dynamics of 
the static and nonequilibrium behavior of these systems 
leading to novel phenomena such as spin textures, spin 
domains, and complex spin mixing dynamics [l|. 

Many of the interesting phenomena in spinor con¬ 
densates hinge on extreme sensitivity to the relative 
strengths of two-body interactions between the internal 
states; those interactions are typically weak for alkali 
atoms (with the exception of ®^Rb [l^, [1^, and possi¬ 
bly also ^^^Cs and ^Li) and as a consequence the spinor 
physics with strong interactions has largely remained un¬ 
explored. Although the usual Feshbach resonance tech¬ 
nique for tuning the interactions cannot immediately 
be applied to spinor condensates, several proposals for 
achieving such control exist [30l - l36l| . which enable possi¬ 
ble ways to study strongly correlated spinor condensates. 
Evidently, as the interactions become resonant, three- 
body contributions need to be considered. In the strongly 
correlated regime, i.e., when the s-wave two-body scat¬ 
tering lengths associated with the hyperfine spin states 
exceed the typical range of interatomic interactions (the 


van der Waals length, r^dw), Efimov physics becomes 
important and several new aspects have to be 

considered. For instance, the existence of an infinity of 
weakly bound Efimov trimers strongly affects the scat¬ 
tering observables at ultracold energies and can poten¬ 
tially impact the spin dynamics in spinor condensates. 
Recent experimental and theoretical work on the three- 
body parameter (dol - f^ has extended the universal prop¬ 
erties of Efimov states allowing for further quantitative 
predictions of few-body phenomena relevant for spinor 
condensates. 

One of the major differences between the usual “spin¬ 
less” Efimov physics and the effect for spinor systems is 
the presence of multiple length scales in the problem. In 
spinor condensates the atomic energy levels are (2/-I- I)- 
fold degenerate (/ is the atomic hyperfine angular mo¬ 
mentum and rrif = —/.../ its azimuthal component), and 
there exist / -I- I rotationally-invariant s-wave scattering 
lengths [H . In Ref. we explored in depth the Efimov 
physics for / = I spinor condensates and showed how the 
presence of the additional internal degrees of freedom can 
lead to multiple families of Efimov states that strongly 
affect the condensate spin dynamics. In the context of 
nuclear physics, where isospin symmetry plays an impor¬ 
tant role, the work of Bulgac and Efimov has also 
demonstrated a much richer structure for Efimov physics 
when the spin degree of freedom was considered. As we 
will see later in this manuscript, a fundamental difference 
of the spin physics in nuclear systems and spinor conden¬ 
sates is that in a gas phase the atoms can be prepared 
in a single spin substate, which is not an eigenstate of 
well-defined total angular momentum. 

Our treatment begins with a summary of the mul¬ 
tichannel generalization of the adiabatic hyperspherical 
representation using contact two-body interactions. We 
then characterize the Efimov physics for spin-1, -2 and 
-3 by determining the topology of the three-body hyper- 


2 


spherical potentials and analysing the possible families 
of Ehmov states present in each case. Special attention 
is given to the spin -2 ®®Rb condensate due to the natu¬ 
rally large values for the atomic scattering lengths found 
in Ref. [13, 113 • Finally, we determine the important 
three-body mean-field parameters, thus building the fun¬ 
damental blocks for future studies of three-body effects 
in spinor condensates. 

II. ADIABATIC HYPERSPHERICAL 
REPRESENTATION FOR SPINOR SYSTEMS 


The mass-scaled Jacobi vectors (written in the “odd- 
man-out” notation) are given in terms of the atomic dis¬ 
tances by 

= (rj - ri)/d, and = d > (4) 

where d = A set 12 of six hyperrangles describe 

all other internal degrees of freedom associated to the 
system’s internal motion. Treating the hyperradius as an 
adiabatic parameter, the three-body wave-function can 
be expressed as 


Our study of few-body physics in spinor condensates 
begins from the multichannel generalization of the zero- 
range Ferr ni p seudopotential for s-wave interactions. In 
Refs. a suitable generalization was proposed to 

be (in a.u.) 


v{r) = -d^(r)—r, ( 1 ) 

m or 

where S^(r} is the three-dimensional Dirac delta- 
function, m is the atomic mass, and A is a scattering 
length operator which incorporates all the important mul¬ 
tichannel structure of the two-body interactions . For 
spinor condensates we assume the scattering length op¬ 
erator in Fq. m can be defined for each pair of atoms 
as 

X! \F2hMF2i,)aF2b{F2hMF2i,\, (2) 

F2hMF^^ 

where 

|F 2 bMF 2 b) = |(/l/ 2 )A 2 bMF 2 b) = 

is the two-body spin functions of total angular mo¬ 
mentum |/i — / 2 I < F 2 h < /i + /2 and projection 
+ ™/ 2 ) expressed in terms of the Clebsch- 
Gordan coefficients. Therefore, the two-body interac¬ 
tion in Fq. o is diagonal in the spin basis {|cr)} = 
We note that, due to bosonic symmetry 
only the symmetric spin states (Fjb even) are allowed 
to interact with rotationally-invariant scattering lengths 
OFab = {“ 0 , 02 , 02 /}- These scattering lengths set im¬ 

portant length scales in the system. Many-body prop¬ 
erties of the system such as the miscibility of spin com¬ 
ponents are sensitive to the relative strength and sign of 
the scattering lengths [I|. Their strength also signifies 
the appearance of universal three-body physics in the 
scattering observables. 

The three-body problem is solved using a multichan¬ 
nel generalization of the adiabatic hyperspherical method 
via a Green’s function approach developed in Ref. (ssl - 
[stII . As usual, the three-body system is characterized by 
a single length scale, the hyperradius R = {p\-\- 


V E 


where F{R) are the hyperradial wave functions, $(i?; 12) 
the channel functions and {|E)} = { 1 ^ 71 / 1 ,^/ 2 , 771 / 3 )} the 
three-body spin functions, formed by the combination of 
all possible product states, chosen to simplify the formu¬ 
lation. The adiabatic channel functions are eigenstates 
of the fixed-i? hyperangular Schrodinger equation 




A2(f2) -f 15/4 

2pR2 


&E' + (F'|I>(i?,f2)|E) 




$^(i?;f2) = [/,(i?)4>^(i?;f2), (6) 


with corresponding eigenvalues given by the adiabatic 
potentials U(R) which determine the hyperradial motion 
and describe (in conjunction with the nonadiabatic cou¬ 
plings) the possible bound states and scattering proper¬ 
ties of the system. In the above equation, /t = mj^/5 is 
the the three-body reduced mass, A is the grand angu¬ 
lar momentum [59l| , V is the sum of all possible pairwise 
interactions [Fq. ([T])], given by 




i<j 




9 Ik) 


dp 




where the superscript (fc) for the “odd-man-out” signifies 
that the interaction is between particles i and j while 
particle k spectates. In the present study the three-body 
energy levels, E^, are degenerate and we set them to 
zero. 

Equation ([5]) is solved using a hyperangular Green’s 
function and the corresponding Lippmann-Schwinger 
equation for each component of the channel func¬ 
tion. 


$E(i?;f2) =-2pR2^ / df2'GEs(^2,I2') 

x4"^,(R,f2')$E,(i?;f^'), (8) 

where 

44(R,f2) = (EKri,)|E'). (9) 
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The components of the hyperangular Green’s function 
satisfy [A^ — (s^ — 4)]GsE(f^, fi') = <5(17, ff'), where s is 
related to the hyperangular eigenvalue U (R) through 


U{R) = 


s{R)^ - 1/4 


( 10 ) 


(k) 

Evaluating the integral over v^-£,, considering only states 
of total orbital angular momentum L = 0, and solving 
the Lippmann-Schwinger equation reduces to determin¬ 
ing values of s for which the determinant of the matrix 


Q = 


2i/2p 


-1 


( 11 ) 


vanishes . In the equation above the matrix ele¬ 

ment for the M matrices are given by 

cot(s7r/2), 


M«,= 


where A^e' is the matrix element for the two-body scat¬ 
tering matrix written in the three-body spin basis, 

= (E|i«|E') 


-A 


(i) 4sin(s7r/6) 

>/3 sin(s7r/2) ’ 


i = 1, 

1 = 2,3, 


( 12 ) 


= «^’2b 

-F2b.M fji, 


(13) 


with the ith particle as spectator to the interaction. The 
P-)- and P_ matrices represent cyclic and anti-cyclic per¬ 
mutations of the the three-body spin basis 


(F+)ee' 

= (S|Fi23|E') 





|Fi23|TO/i,m/3, 

.m%) 




(14) 

(F-)es' 

= (E|Fia2|E') 





|Fi32|TOj^,m/^, 



,rrj 


(15) 


The spin basis {|E)} was chosen to to give simple ex¬ 
pressions for the components of Q which were written 
in the odd-man-out notation. In this form, however, Q 
does not account for any symmetry property in the sys¬ 
tem. This is fixed by the unitary transformation SQS'^ 
with 


(^)ee' = (S| ^ \F3bMF,,iF2b)){F3^Mp,,{F2b)\ |S'), 

VF2b / 

(16) 

where the three-body spin functions of total angular mo¬ 
mentum |P 2 b — /I < Tab < \F 2 h + /I and projection 
M= MF^-b F ITT-f are given by 

\F3hMF3b{F2h)) = |(/l/2/3)-F3bAf_F3t(F2b)) = 

{F2hMF2b fsmfs I Fab Mf^i, ) \ Fab) I /a w/ 3 ). 

( 17 ) 


Therefore, evaluating the det[F(5F^]=0 results in a tran¬ 
scendental equation whose roots s{R) determine the 
channel functions $(P; fl) and the three-body potentials 
U{R) from Eq. (flUl) with well define hyperfine angular 
momentum (Fab and MF^b) and permutation symmetry. 
In the usual spinless problem, the solution for three iden¬ 
tical bosons in the limit P/a —>■ 0 (a been the scattering 
length) gives a lone imaginary root sq ~ 1.0062i, which 
produces an attractive three-body potential when sub¬ 
stituted in Eq. OT. This potential supports an infinite 
number of bound trimers characteristic of the Efimov ef¬ 
fect. In the spinor case elaborated above, the degeneracy 
of the hyperfine spin-manifold yields fundamentally dif¬ 
ferent three-body physics as demonstrated in the follow¬ 
ing two sections. 


III. SPIN-1 SYSTEMS 


In this section, the results for spin-1 from Ref. are 
summarized, focusing on the structure of the three-body 
potentials and the families of Efimov states. The impact 
of Efimov physics on the three-body scattering observ¬ 
ables and the many-body physics is a subject which is 
delayed until Sec. V. In order to explore the solutions s 
of det[P(5<S^]=0 and the physics resulting from them, we 
assume the relevant scattering lengths for the problem, 
OFab = {ao) 02 ) < 22 /}, to differ greatly in magnitude 

from each other. By doing so, the values for s with a 
range, say, |oq,| <C P <C lo/jj is constant and the corre¬ 
sponding values can be easily listed and emphasize the 
attractive (s imaginary) or repulsive (s real) 1/P^ inter¬ 
action that characterizes the Efimov physics. When this 
condition of highly different scattering length magnitudes 
is not satisfied, one must in general return to the exact 
solutions of the above transcendental equation. In the 
regimes considered here, there will exist typically various 
regions in which a detailed analysis must be performed 
(see Fig. [T|). Here this analysis is carried out for every 
possible case. 

Following our analysis in the previous section we can 
derive the transcendental equations that determine s^, 
and consequently the three-body potentials Up{R) [see 
Eq. ([Toll ]. For / = 1 atoms, the possible values for Fsb 
are 0, 1, 2 and 3. An analysis of the three-body spin 
functions (HZl) for Fab = 0 [^, however, shows that such 
states are fully anti-symmetric and are non-interacting. 
[The potential model in Eq. ([T]) only captures s-wave 
interactions.] For Fab = 1, however, the spin functions 
are fully symmetric and the values for s are determined 
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FIG. 1: Schematic representation of the relevant regions in 
R of the three-body potentials U{R). The shaded area, R < 
f'vdw, represents the region where the three-body potentials 
are expected to be non-universal. For all other regions, the 
potentials are universally given by the form in Eq. (iia with 
s{R) = s. 


by solving the transcendental equation. 


3* (oq -b 02)5 cot(§s) 32 ( 0002 ) 3 ^ cot(5s)^ 


2^R 

25 (2oo -b 02 ) sin(|s) 
33 sin(|s)i? 

16 aoa2sin(^s)^ ^ 

35sin(fs)2i?2 “ 


2i?2 

2 00025 cot(^s) sin(^s) 

sin(^s)i?2 

(18) 


(Imaginary roots can be obtained by mapping s —>■ is.) 
The above transcendental equation depends on both two- 
body scattering lengths, oq and 02 . For Fab = 2 and 3 
the spin states are symmetric and mixed-symmetry states 
[ 5 ^, respectively, and the interaction strength s is ob¬ 
tained, respectively, through, 

3^/^O2SC0t(|s) ^ 2 2 ^/ 2 q 2 s sin(|s) 

2 i/ 2 i? + 31/4 sin(f s)i? 

and 

3^/'^02SC0t(|s) ^ 4 2^/2 q2s sin(|s) 

2 i/ 2 i? + 31/4 sin(f s)i? 


(19) 


( 20 ) 


Note that for Fab = 2 and 3 the transcendental equa¬ 
tions depend only on 02 . Table U lists the solutions of 
Eqs. (I18II - ((^D11 for the regions in R in which the value of 
s is constant. Note that we also list the relevant scat¬ 
tering lengths for each value of Fab. For the cases where 
an imaginary root exists, i.e., when Efimov states are al¬ 
lowed, we list both the imaginary root and the lowest 
real root. For cases in which an imaginary root does 
not exist, we list only the lowest real root. It is inter¬ 
esting noting that when R is much smaller than all rele¬ 
vant scattering lengths, R <C |a{o, 2 }|j the imaginary root. 
So ~ 1.00621, is the same than the usual three identical 
boson problem. For the regions |ao| R |a 2 | and 


|a 2 | <C F <C |ao|, i.e., when the oq and 02 are not effec¬ 
tively resonant, respectively, the possible roots can be dif¬ 
ferent than the ones obtained for the usual bosonic prob¬ 
lem. These new roots result from the fact that the spin 
functions [Eq. (ITTI) ] are superpositions of single particle 
product states (as determined by the angular momentum 
addition algebra). It is also interesting noting that for all 
cases we have explored, mixed-symmetry states (Fab = 2 
for / = 1) display a stronger repulsive barrier whenever 
R is larger than all scattering lengths (s = 4) in compar¬ 
ison to symmetric states (s = 2). This implies that some 
three-body scattering observables should be suppressed 
for mixed-symmetry states. This will, in fact, simplify 
some of our analysis in Section V. 

Figure [2] shows the three-body potentials for / = 1 
as a result of solving Eq. m for the allowed values of 
Fab. The scattering lengths chosen illustrate only one 
of a pair of possibilities for the relative strengths of the 
scattering lengths. The hyperradius is written in terms 
of the typical range of interatomic interactions, rvdw. 
The totally symmetric Fab = 1,3 states, which are the 
only states not suppressed by centrifugal barriers, are 
associated with a pair of attractive potentials for R <C 
{ 0 - 0 , 02 }. This pair is associated with sq ~ 1.00621, and 
allow for the coexistence of two families of Efimov states 
(represented in Fig. [2] by the horizontal solid and dash- 
dotted lines.) As the hyperradius exceeds oq, the Fab = 
1 attractive potential becomes an atom-dimer channel 
for collisions between the dimer |F 2 b = = 0) 

with energy — l/mog and a free atom in the jm/ = 0) 
state. In the intermediate region oq F <C 02 , only 
the Efimov potential in the Fab = 3 state remains and is 
associated with sg ~ 1.00621. For 02 <C F, the Fab = 3 
Efimov potential joins repulsive barriers from Fab = 1,2 
in an atom-dimer channel |F 2 b = 2,m/ = —1,0, l)|w/ = 
— 1,0,1), which introduces the interesting possibility of 
studying atom-dimer spin mixtures [s^ j. 

As a final note, if a small (but finite) field is applied, 
the Efimov states in Fab = 1 and 3 can interact creat¬ 
ing an overlapping series of states. Such controllability 
can produce ultra-long-lived states whose presence 
can affect the spin dynamics of the condensate depend¬ 
ing on the short-range physics . In the case of higher 
spin analyzed in the following section, many families of 
Efimov states with additional novel roots appear, along 
with the presence of overlapping series of Efimov states 
in the absence of an applied field. 


IV. SPIN-2 AND -3 SYSTEMS 

In this section we go beyond the analysis from Ref. [1^ 
and detail the three-body potentials and Efimov states 
for / = 2 and 3 spinor condensates. For spin-2 there are 
three rotationally invariant scattering lengths (ag, 02 , a 4 } 
and Fgb = 0...6. The three-body spin states [Eq. (fT71) ] for 
Fgb = 0, 2, 3, 4, and 6 are symmetric while for F^h = 1 
and 5 they are of mixed-symmetry. Shown in Table I are 
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R (lO^rvdw) 


(b) (R=ao) {R=a 2 ) 



FIG. 2: Fab = 1 (red solid line), 2 (green dashed line), and 3 (blue dash-dotted line) hyperspherical adiabatic potentials for 
/ = 1 atoms with ao = lO^rvdw and 02 = lO^rvdw- (a) For R < { 00 , 02 } (shaded region) two attractive potentials exist 
(both with So ~ 1.0062i), allowing for two families of Efimov states, and for R > ao, one of these potentials turns into an 
atom-dimer channel |F 2 b = 0, Mfj,, = 0)|/ = l,mf = 0). (b) For oo < i? < 02 (shaded region), only one family of Efimov 
states exists (so « 1.0062i), and for i? S> 02 three (asymptotically degenerate) potentials describe atom-dimer channels, 
|F2b = 2, Mf,, = -1, 0,1)1/ = 1, m/ = 1, 0, -1). 


the relevant roots Sjy and relevant scattering lengths for 
each Fsb state. Similar to the / = 1 case, for / = 2 
the roots for the regions in which i? is smaller than all 
relevant scattering lengths correspond to the roots for the 
usual identical bosons problem. All other regions display 
new roots. For instance, within the region jodj <C i? <C 
|a{o. 2 }| novel imaginary roots sq ~ 0.3788f (Fab = 2) and 
So ~ 0.5528f (Fab = 3) can be found. Analogously to the 
discussion of the spin configurations in Section III, the 
presence of these roots can be explained by recognizing 
that the three-body spin functions are not simple product 
states. 

Although typically the scattering lengths for most al¬ 
kali species are small (comparable to rvdw)j one no¬ 
table exception is ®®Rb with / = 2. For that system 
an ~ —8.97rvdw, 02 ~ —6.91rvdw, and 04 Ri —4.73rvdw 
1^ . Recent experimental and theoretical advances 
with alkali atoms have determined that the value of the 
scattering length at which the first Efimov state appears 
is approximately — lOrvdw [dOl - ls^ . Therefore, for / = 2 
®®Rb spinor condensates, it is not unreasonable to expect 
that the Efimov physics can strongly impact the spin dy¬ 
namics and many-body physics. In Fig. |31 the hyper¬ 
spherical adiabatic three-body potentials for spin -2 ®®Rb 
are shown. As one can see, the potentials for Fab = 0, 2 , 
3, 4 and 6 states support an attractive Efimov potentials 
for distances smaller than and comparable to the rele¬ 
vant scattering lengths, but forming a repulsive barrier 
as the hyperradius exceeds those distances. The poten¬ 
tials for the mixed-symmetry states. Fab = 1 and 5, are 
completely repulsive and support no Efimov states 

Finally, for three identical spin-3 atoms there are four 
rotationally invariant scattering lengths { 00 , 02 , 04 , 06 } 
and Fab = 0, ...,9. For spin-3, the Fab = 0 spin state 
is fully anti-symmetric, while Fab = 2 and 8 states are 



FIG. 3: Fsb = 0 (red solid lines), 1 (red dot-dashed lines), 2 
(green dashed lines), 3 (blue dashed lines), 4 (purple dashed 
lines), 5 (cyan dot-dashed lines) and 6 (black dot-dashed lines) 
hyperspherical adiabatic potentials for / = 2 ®®Rb atoms with 
ao ~ —8.97rvdw, 02 ~ —6.91rvdw, and 04 Ri —4.73rvdw. 
Note that the potentials for the mixed-symmetry states, 
Fsb = 1 and 5, do not display an Efimov potential. 


mixed-symmetry states and all others being fully sym¬ 
metric states. Candidates for interesting Efimov physics 
(i.e., some of the relevant scattering lengths are large) are 
®®Rb, ^^^Cs and ^Li. In Table I, the values of s^, for all 
values of Fab and all hyperradial ranges are presented for 
the first time. Where the spin-2 analysis produced only 
a few novel imaginary values of Sj/, there are many more 
families Efimov states for spin-3. Of particular note is 
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TABLE I: Values of s^, relevant for /=!, 2 and 3 spinor condensates covering all possible regions of R and for different magnitudes 
of the relevant scattering lengths. We list the lowest few values of for each Tab and their multiplicity (superscript) if greater 
than one. 


(/ = 1) 

Relevant aF2b 

^3b — 1 

{ao, 02 } 

^3b — 2 
{«2} 

Fsh = 3 
{«2} 





R‘C|a,{o, 2 } 1 

1.0062i, 2.1662 

2.1662 

1.00622,4.4653 





|aol<Sfl<g|a 2 | 

0.7429 

2.1662 

1.00622,4.4653 





io 2 i<gJi<giaoi 

0.4097 

4 

2 





i?;S>|a{o, 2 } 1 

2 

4 

2 





(/ = 2) 

^3b = 0 

Fsb = 1 

Fsh — 2 

f^3b — 3 

Fsh = 4 

-fab — 5 

Fab = 6 

Relevant aF2b 

{^2} 

{“ 2 } 

{ao,a 2 , a4} 

{a 2 , 

{a 2 ,a4} 

{ 04 } 

{ 04 } 

R<^\Q'{0,2,4} I 

1.00622, 4.4653 

2.1662 

1.0062i, 2.1662^^) 

1.00622, 2.1662 

1.00622, 2.1662 

2.1662 

1.00622,4.4653 

|oo|<ii<|a{2,4}l 

1.00622, 4.4653 

2.1662 

0.49050 

1.00622, 2.1662 

1.00622, 2.1662 

2.1662 

1.00622,4.4653 

io2j<gii<gja{0,4} 1 

2 

4 

0.74732, 2.1662 

1.1044 

0.66080 

2.1662 

1.00622, 4.4653 

io4j<ii<ja{(j_2}i 

1.00622, 4.4653 

2.1662 

0.37882, 2.1662 

0.55282, 3.5151 

0.52186 

4 

2 

k{ 0 , 2 } l<SJJ<|o4i 

2 

4 

0.97895 

1.1044 

0.66080 

2.1662 

1.00622, 4.4653 

i“{0,4} i<Sii<gjo2j 

1.00622, 4.4653 

2.1662 

1.3173 

0.55282, 3.5151 

0.52186 

4 

2 

k{2,4}i<Ji<iooi 

2 

4 

0.68609 

2 

2 

4 

2 

-R^|Q'{0,2,4} I 

2 

4 

2 

2 

2 

4 

2 

(/ = 3) 

^3b — 1 

^3b — 2 

7^3b — 3 

-fab — 4 

Fab = 5 



Relevant aF2b 

{a2,04} 

{02,04} 

{ao, 02 , 04 , ae} 

{a2,a4, ae} 

{a. 2 , 04 , ae} 



-R<C|a,{0,2,4,6} 1 

1.00622, 2.1662 

2.1662(^) 

1.0062h^L 2.1662'^) 

1.0062i, 2.1662*^) 

1.0062i, 2.1662(^> 



|oo|<ii<|a(-2,4,6}l 

1.00622, 2.1662 

2.1662(^' 

1.00622, 0.3420 

1.0062i, 2.1662<^) 

1.0062Z, 2.1662<^> 



|o2|<gii<g|a{o,4,6} 1 

0.6608 

2.1662 

1.00622, 0.5469 

0.8754 

0.2588 



ja4j<CR<cja-.{-0,2,6} 1 

0.5219 

2.1662 

1.00622, 1.1901 

0.93522, 2.1662 

0.66782, 2.1662 



|a6|<ii<|a{0_2,4}l 

1.00622, 2.1662 

2.1662(^> 

1.00622, 0.41122 

0.43092, 2.1662 

0.3351i 



I‘3'{o,2} 1 

0.6608 

2.1662 

0.65212, 1.0098 

0.8754 

0.2588 



I‘3'{o,4} 1 

0.5219 

2.1662 

0.6080i, 1.6045 

0.93522, 2.1662 

0.66782, 2.1662 



l“{ 0 , 6 } l<^i<|a{2,4} 1 

1.00622, 2.1662 

2.1662 

0.9666i, 1.1343 

0.43092, 2.1662 

0.33512, 2.1662 



l“{2,4} l<SJi<|a{0,6} 1 

2 

4(2) 

0.5858i. 1.7691 

1.0111 

0.9552 



l“{ 2 , 6 } i<Ji<sia{0,4} 1 

0.6608 

2.1662 

0.9100i, 1.1558 

1.6693 

1.2189 



io{4_6} i<gii<sia{0,2} 1 

0.5219 

2.1662 

0.42892, 1.2015 

0.0803i, 3.3774 

0.8199 



k{0,2,4} l*CR^|n 6 1 

2 

4 

0.9984 

1.0111 

0.9552 



1 <^{ 0 , 2 , 6 } |■^-R■^|a4 1 

0.6608 

2.1662 

0.64152, 2 

1.6693 

1.2189 



|o{0,4,6} 1 

0.5219 

2.1662 

0.6392 

0.0803i, 3.3774 

0.8199 



|a-{2,4,6} |*C-R*c|ao 1 

2 

4 

0.7819 

2 

2 



-R!!^|a,{o,2,4,6} 1 

2 

4(2) 

2(2) 

2 

2 



(/ = 3) 

^3b = 6 

Fab = 7 

-fab — 8 

f^3b — 9 




Relevant aF 2 b 

{04, ae} 

{04, ae} 

{“el 

{“el 




-^*^ 10 ,( 0 , 2 ,4, 6 } 1 

1.00622, 2.1662 

1.0062i, 2.1662 

2.1162 

1.00622,4.4653 




|ao|<ii<|a(;2,4,6}l 

1.00622, 2.1662 

1.0062i, 2.1662 

2.1162 

1.00622,4.4653 




|a2|<gfl<g|a{0,4,6} 1 

1.00622, 2.1662 

1.00622, 2.1662 

2.1162 

1.00622, 4.4653 




|a4|</J<|a{(],2,6} 1 

1.1329 

0.6372 

2.1162 

1.00622, 4.4653 




ia6j<fl<ja{o,2,4}i 

0.58422, 3.5329 

0.5491 

4 

2 




k{ 0 , 2 } l<^i<S|a{4,6} 1 

1.00622, 2.1662 

1.00622, 2.1662 

2.1162 

1.00622, 4.4653 




i“{0,4} i<^i<ia{2,6} 1 

1.1329 

0.6372 

2.1162 

1.00622, 4.4653 




1 O'{ 0 , 6 } i‘=^-^‘=^io{2,4} 1 

0.58422, 3.5329 

0.5491 

4 

2 




l“{2,4} l<^i<|a{0,6} 1 

1.1329 

0.6372 

2.1162 

1.00622, 4.4653 




l“{ 2 , 6 } i<JJ<ia{0,4} 1 

0.58422, 3.5329 

0.5491 

4 

2 




ia{4^6} i<^i<ia{ 0 , 2 } 1 

2 

2 

4 

2 




k{0,2,4}I^^J^|a6| 

1.1329 

.6372 

2.1162 

1.00622, 4.4653 




|O{0,2,6 } |*^-R*^|o4 I 

0.58422, 3.5329 

0.5491 

4 

2 




10{0,4,6} 1^2 1 

2 

2 

4 

2 




|o{2,4,6} l*CR^|ao 1 

2 

2 

4 

2 




-f^^|o{0,2,4,6} 1 

2 

2 

4 

2 





the appearance of overlapping Efimov series for Fab = 3 
when R |a{o, 2 . 4 , 6 }l and |a{6} <C i? <C |a{o, 2 , 4 }|- Even 
at zero applied field these overlapping states have the 
potential of being ultra-long lived depending on the 
short-range physics [s^. Analysis of these overlapping 
states is beyond the scope of the present work. The fi¬ 
nal section concludes with an analysis of the impact of 
Efimov physics on the three-body observables and many- 
body physics. 


V. TWO- AND THREE-BODY MEAN-FIELD 
PARAMETERS 

In this section we discuss the important two- and three- 
body mean-field parameters for spinor condensates and 
their connection to Efimov physics. It is important to 
notice that in typical spinor condensate experiments the 
system is prepared at some small but finite value of the 
magnetic field. In that case, the Zeeman interaction 
breaks the degeneracy of the atomic levels and atoms 
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can be prraared in a single hyperfine substate, typically 
w/ = 0 [l|. After this preparation state, the magnetic 
field is quickly ramped off making now all atomic lev¬ 
els degenerate and, consequently, making other hyperfine 
substates energetically available. As a result, the system 
evolves from a state where all atoms are in the m/ = 0 
to a more complicated mixture of spins. Evidently, this 
spin dynamics is primarily controlled by the interatomic 
interactions since the initial state is not an eigenstate of 
the total hyperfine angular momentum. 

Therefore, in order to connect with many-body treat¬ 
ments and identify the important interaction terms that 
control the spin dynamics it is necessary to rewrite the 
rotational invariant interactions in a suitable form. One 
such way was elaborated in Refs. [Slil for two-body 
interactions. In this case, the two-body scattering length 
operator A [see Eq. ©] can be written as 

^ ( 21 ) 

^2b 


where PF^b = lf' 2 hMF^J(F 2 hMr^^^l are projec¬ 

tion operators into a two-body total hyperfine state. In 
Ref. [ig it was shown that the relation, 


(/■A,y=E 

^2h 


P2b (P2b + 1) 


-/(/ + 1 ) 




( 22 ) 


allows us to rewrite A as 


/ 

^ ^ “2b 

n—0 


(n) 


{h-h)\ 


(23) 


where the Q! 2 b’s are linear combinations of the scattering 
lengths {ao,a 2 , ..., 02 /}. 

The above formulation allows us to rewrite the scat¬ 
tering length operator A in Eq. (l23t for / = 1 atoms in 
terms of two parameters. 


A = o™ -t 


(/i • fV) 


where 


(0) OO + 202 j (1) 02 — Oq 

‘^2b =-5- “2b =-5-’ 


(24) 

(25) 


representing a direct interaction terms and an spin- 
exchange term, respectively. It has been shown in Refs. 
[il[i3 that these two parameters do characterize impor¬ 
tant phases of the gas. For instance, depending on the 
sign of , the ground state of the spinor condensate is 
antiferromagnetic ( 02 b ^ ^) ferromagnetic (o^b ^ *^) 
For / = 2 , the scattering length operator A 
in Eq. (|23p can be rewritten in terms of three parameters 


a 

a 

a 


( 0 ) 

2ao 

802 

2 b 

~ 5~ 


( 1 ) 

_ ao 

202 

2 b 

30 ’ 


( 2 ) 

_ ao 

02 0 

2 b 

~3d~ 

21 7 


a4 

70 


9a4 

9o4 

TiT’ 


(26) 


We note that in the literature the A operator for / = 2 
has been rewritten in terms of the projector Pq allow¬ 
ing for a structure similar to d-wave BCS superfluids 
[6^ . For / = 2, there are three categories of ground 
state phases depending on the relative strengths of the 
scattering lengths: cyclic, polar, and ferromagnetic (see 
Ref. [i^.) The cyclic phase arises from the additional Po 
term which describes scattering of a pair of atoms into 
the F 2 b = 0 singlet state introducing an additional order 
parameter describing the formation of singlet pairs in the 
system. For / = 3, the four parameters that characterize 
A [Eq. (|2^ ] are: = 9ao/35 — 4a2/7 -I- 48604/385 -I- 

406/77, = 9oo/70 - 17o2/63 -f 8 I 04 / 77 O -f 25o6/693, 

= — 00 / 3 I 5 -I- 502/378 — 604/385 -I- 23o6/4158, and 

= —oo/630 -I- 02/378 — 04 / 770 -I- 06/4158. A was 
rewritten in Ref. [63| including Pg and rewritten in 
Ref. [6^ including a nematic tensor for the dipolar spinor 
condensate ®^Cr using the terminology of liquid crystals 
for the many ground state phases of the system at zero 
and nonzero fields which we do not go into further 
There are, however, many ways to rewrite A in the 
literature that better suggest treatments similar, for in¬ 
stance, to magnetism, BCS theory, and liquid crystal the¬ 
ory. For the three-body mean-field contributions, how¬ 
ever, we choose to simply extend the above analysis and 
in an analogous form to Eq. (l23)l . in light of the lack 
of suggestive three-body treatments from other theories 
to our knowledge. Therefore, we define the three-body 
scattering length operator as 

^b = E«3b^‘’^^^’3b> (27) 

^3b 


where 

7’F3b= E |7"3bMF3,(E2b))(F3bMF3,(E2b)|. (28) 

MF3b Fab 


In the equation above, asb is the usual three-body 
scattering length (units of length"^) as defined in Refs. 
[^ [6416^ for the allowed values of Fsb. Now, using the 
usual relations for angular momentum addition and the 
orthogonality of the projection operators, one can obtain 
the following relation. 


Ea* =E 

. i<3 ) Fsh 


FsAFsh + l) 3/(/ + l) 


Ff,, 


(29) 

which allows us to write the three-body scattering length 
operator Agb as 


Af3b-1 / 

^3b = E “ 3 ^ 

n=o y 

where Aab is the number of relevant three-body scat¬ 
tering lengths, i.e., the three-body scattering length for 




i<j 


(30) 











the values of whose spin functions are fully symmet¬ 
ric. (Anti-symmetric and mixed-symmetry states corre¬ 
spond to higher partial waves and are not considered in 
the present work.) Therefore, the three-body mean-field 
contributions to the spin dynamics is now determined in 
terms of the parameters osb, which are linear combina¬ 
tion of the physical three-body scattering lengths for all 
allowed values of Fab. 

For / = 1 atoms, only the three-body scattering 
lengths for Fab = 1 and 3 contribute to the interaction 
(see Table I), allowing us to write the three-body scat¬ 
tering length operator Aab in Eq. (15(1)) in terms of only 
two parameters. 


^3b = ttab + ttab ’ 


. i<j 


where 


(31) 


Q (1) , O (3) (3) (1) 

(0) _ 3a3j^ -I- (i) _ 

'^3b — 5 g , 


(32) 


representing a direct interaction term and an spin- 
exchange term, respectively. This is in close analogy to 
the form of the two-body scattering length operator for 
/ = 1 atoms [see Eq. (IMl) ]. For / = 2 atoms, only the 
states with Fab = 0, 2, 3, 4, and 6 contribute in the inter¬ 
action, resulting in the following parameters for Eq. ()30)) , 


v(o) 


2 a 


(0) 

3b 


2a 


( 2 ) 

3b 


35 


+ 


3a, 


(3) 

3b 


(0) 

( 1 ) ^ 

30 


23a 


( 2 ) 

3b 


29a: 


+ 

(3) 

3b 


243a, 


(4) 

3b 


385 

513a 


(6) 

Qsb 

385’ 


(4) 

3b 


126 


60 


v( 2 ) 


(3) 

= - 


(4) _ « 
“3b 3780 


.( 0 ) 


.( 2 ) 


,(3) 



1 

^"^^3b 1 

1260 

' 126 


360 

“3b 1 

( 2 ) 

“3b 1 

(3) 

'^^3b 

945 

1134 

540 

1540 

,( 0 ) 

'3b 

“3b 

^3b 

^3b 


1540 

117a 


“3b 

990’ 


(4) 

3b 


3080 


+ 


( 6 ) 

3b 

770’ 


a. 


17a, 


(6) 

3b 


62370’ 


a 


(6) 

3b 


1134 1080 3080 62370 


(33) 


For / = 3 atoms, the degree of complexity increases 
rapidly where now the states with Fab = 1, 3, 4, 5, 6 , 
7, and 9 contribute in the interaction (see Table 1) and 
analog expressions for oab can be obtained. 

In order to estimate the two- and three-body mean- 
field contributions and their relative importance, we an¬ 
alyze the two- and three-body coupling constants 


(0 ^TT (i) 

52b = -« 2 b 


and 


(i) _ 31 / 2 !^ (i) 

53b — *4 ^ «3b > 


(34) 


respectively [ 6 ^. This way, we can estimate the two- and 
three-body mean-field energies simply as ng 2 h and n^gzh, 
respectively, with n being the atomic density. This makes 
clear that relative importance of three-body effects will 
depend on the density (the denser the gas the more im¬ 
portant three-body contributions are.) For instance, as 


pointed out in Ref. [H^, the / = 1 ferromagnetic and 
anti-ferromagnetic phases can be affected by three-body 
physics whenever the three-body mean-field energy ex¬ 
ceeds the two-body mean-field energy, i.e., 

«^l53b I > «l52b I, (35) 


and are of opposite sign. It is interesting noting that 
for most of the alkali species the two-body scattering 
lengths are typically small, implying that the two-body 
direct and spin-exchange mean-field energies are also 
small opening up ways to explore the importance of 
three-body effects for the system. The fact that the two- 
body scattering lengths are small does not necessarily 
imply that the three-body scattering lengths, a 3 b, are 
also small. A notable case is the / = 1 ®^Rb spinor 
condensate where not only are the scattering lengths 
small, but they are also approximately the same. For 

this case, ao ~ 1.23rvdw and a 2 = 1 . 21 rvdw [13, 

implying an extremely small spin-exchange energy term 
~ —5.6 X 10“®rvdw [see Eq. ([25])], and the three- 
body spin-exchange can be important. Determining the 
three-body contributions (the actual value of the three- 
body scattering lengths) for such cases is extremely chal¬ 
lenging since it would require a full numerical calculation 
including realistic two- and three-body interactions. 

In the strongly interacting regime, not only do the 
three-body correlations become important, but they also 
become universal [5^. The Efimov physics for the spin 
problem displays many features that can allow for an in¬ 
dependent control of both two- and three-body physics 
opening up ways to strongly modify the spin dynamics 
in the condensate. From the analysis above, it is clear 
that the impact of Efimov physics to the spin dynamics 
is made through its role in three-body elastic processes. 
In general, Efimov states are manifested in the three- 
body scattering observables through log-periodic inter¬ 
ference and resonant effects. In fact, for / = I atoms 
(see Ref. Hi), a simple WKB model can be used to 
determine the general form of the three-body scattering 
length. For (/ = 1), an interesting case emerges when 
both ao and a 2 assume large and negative values. In that 
case, if |ao| |a 2 | an analytical expression for the two 
relevant three-body scattering lengths (Fab = 1 and 3) 
can be determined and are given by 


“3b “ 


“3b ~ 


a — j3 tan 

fsoln^) 


V “i /. 

a — P tan 

^soln^) 


k “3 J \ 


f ^ 

\ao 


0.82 

4 I 4 
“O + 7“0: 


a 


4 

2 ; 


(36) 

(37) 


where sq ~ 1.0062, and a, /3 and 7 are universal con¬ 
stants that can be determined from numerical calcula¬ 
tions for each value of Fab. [The factor (a 2 /ao) above 
originates from the regions where the three-body poten¬ 
tials are repulsive —see Table I.] For negative values of 
ao and a 2 , the Efimov physics is manifested in elastic 
processes when a Efimov state becomes bound. In the 
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above equations for < 0 are the values 02 < 0 
where a Fsb = 1 and 3 Efimov resonance occurs. Near 
such values of , the three-body scattering length osb 
can be extremely large, assuming both positive and neg¬ 
ative values, allowing for control of the spin dynamics 
via the spin-exchange term in Eq. (1321) . In reality, since 
three-body losses are present and Efimov states have a 
finite lifetime, osb does not strictly diverge and, in fact, 
is a complex quantity whose real and imaginary parts 
describe elastic and inelastic collisions—one can intro¬ 
duce loss effects in the expressions above (and below) by 
adding an imaginary phase term irj in the argument of 
the tangent, where ri is the so-called three-body inelas- 
ticity parameter sl, Nevertheless, for typical values 
of 77 [^, [33 one can still expect substantial tunability of 
the three-body spin-exchange dynamics. 


For / = 2 atoms, evidently, the complexity in de¬ 
termining the three-body spin-dynamics increases. As 
we mentioned above, an interesting case is presented for 
/ = 2 ®®Rb atoms, where all the relevant two-body scat¬ 
tering lengths are large and negative (oq ~ —8.97rvdWj 

02 ~ —6.91rvdW: and 04 ~ —4.73rvdw [13, H^) and 

close to the value in which identical bosons display 
an Efimov resonance, — lOcvdw [dOl - l^ . This implies 
that Efimov resonances might enhance three-body spin- 
exchange. Since the magnitude of the scattering lengths 
in this case are comparable, it makes it difficult to obtain 
an accurate universal expression for the three-body scat¬ 
tering lengths. Nevertheless, it is instructive to analyze 
the results obtained by assuming |ao| |a 2 | lodj. In 
this case, the relevant three-body scattering lengths are 
given by 


(0) 

^3h = 


( 2 ) 

= 


(3) 

%h = 

(4) 

= 

(6) 

^Sb = 


a — P tan ( Sq In ^ 


a — P tan ( Sq In -|- Sg In — 
02 0,4 

a — P tan ( Sg In -|- Sg In — 
0,0 0^4 


\ 1.37 

“2 \ 4 , 4 

— Og -f 7O0 

0-0 J 


. 614 

a — P tan ( Sg In ^ 
04 

, CL 4 

a — P tan I sg In ^ 


X 1.04 

®4 \ 4 4 

— 02+ 702 

02 J 


4 

04 


(38) 

(39) 

(40) 

(41) 

(42) 


where so ~ 1.0062, Sg ~ 0.3788, and Op^^ < 0 are the val¬ 
ues of the two-body scattering lengths where a F 30 Efi¬ 
mov resonance occur. Note that now two distinct families 
of Efimov states (associated to sg and Sg) can affect the 
spin dynamics. Although we expect the values for Op^^ 
to be close to -lOrydw an analysis similar to the one for 
identical bosons |50l | is necessary to precisely determine 
such values. 

For higher spin the three-body contribution introduces 
additional complexity on the manifold of possible ground 
state phases of the spinor condensate. The full study 
of the impact of three-body contributions on the many- 
body phases of spinor condensates remains however the 
subject of future investigation. 


the spin dynamics of the condensate. For spinor sys¬ 
tems multiple families of Efimov states can exist, lead¬ 
ing to novel scaling laws and interesting spin dynamics. 
For small applied field, the families of Efimov states of 
different spin states can couple and may produce ultra 
long-lived states. Curiously, in the case of spin-3 atoms 
there exist overlapping Efimov series even at zero applied 
field. The three-body hyperspherical adiabatic potentials 
for ®®Rb are analyzed to illustrate possible strong effects 
due to Efimov states. Finally, the existence of Efimov 
states can impact the three-body scattering observables 
and many-body phases of the system, and offer novel 
regimes for the spinor system. 


VI. CONCLUSION 
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